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Abstract 

A plasma blob is modeled as consisting of two homogeneous spheres of equal radius and equal 

^I^ ' but opposite charge densities that can move relative to each other. Relative translational and 

Ph. 

rotational motion are considered separately. Magnetic effects from the current density caused by 

^ ■ the relative motion are included. Magnetic interaction is seen to cause an inductive inertia. In the 

Qh, relative translation case the Coulomb attraction, approximately a linear force for small amplitudes, 

Q ' causes an oscillation. For a large number of particles the corresponding oscillation frequency will 

C/3 ■ 

^~>' not be the Langmuir plasma frequency, because of the large inductive inertia. For rotation an 

^. 

Qh' external magnetic field is included and the energy and diamagnetism of the plasma in the model 
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is calculated. Finally it is noted how the neglect of resistivity is motivated by the results. 



I. INTRODUCTION 

In order to study qualitative effects of magnetism in plasma dynamics a very simple 
model is introduced. Two overlapping homogeneous spheres of equal radii, and of equal but 
oppositely signed charge densities, are assumed to move, relative to each other, with negli- 
gible dissipation (resistivity) under the influence of electric and magnetic interaction. The 
neglect of dissipation is motivated at the end. We first treat relative translation (oscillation) 
and then relative rotation. 

Many years ago Tonks and Langmuir [l| carefully derived an equation of motion for a 
collective translational motion of the electrons relative to the positive ions. Their derivation 
seems to indicate that there should be a universal frequency for this mode. 
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^p = \hr^ (1) 



m 



the plasma frequency, depending only on the electron number density n. If a large number 
of electrons move relative to the positive ions one gets a large current and thus it seems 



I2 



as if magnetic effects should affect the result. Bohm and Pines |2| studied the influence of 
magnetic interaction on plasma modes but they did not come up with any explicit correction 
to Up. In the textbook by Goldston and Rutherford ^ the absence of magnetic effects is said 
to be due to a displacement current that compensates for the electron current. In this article 
the problem is approached from a very fundamental starting point: the relevant Lagrangian 
density. The conclusion is that the frequency is lowered by the large inductive inertia [see 
Eq. mi 

We then study the relative rotation of the two charged spheres. If magnetic interaction 
is neglected the kinetic energy is simply determined by angular momentum and moment 
of inertia. When magnetic interaction is included the kinetic energy for a given angular 
momentum is much smaller. The reason for this is again that the effective moment of 
inertia will be dominated by inductive inertia. By adding an external magnetic field to the 
model we can calculate the response of our model plasma and it turns out to be diamagnetic. 



II. SEPARATION OF OVERALL TRANSLATION AND ROTATION 



The kinetic energy of any system of particles 



can be written 



T = E \^^^l (2) 



T = -Mv"^ + -nm + r (3) 

2 2 ^ ' 



where M is total mass, J is the (instantaneous) inertia tensor, V center of mass velocity, 
and fi is a well defined average angular velocity j^, |5|. V is chosen so that p = MV is 
the total momentum of the system and ft so that X = JJl is the total (center of mass) 
angular momentum. T' is the kinetic energy of the particles relative to the system that 
moves with the center of mass velocity and rotates with the average angular velocity. We 
call this the co-moving system. One can introduce generalized coordinates so that there are 
six degrees freedom describing center of mass position and of average angular orientation, 
while T' depends on the remaining 3A^ — 6 generalized coordinates. 

Now consider a blob of plasma that consists of two spheres of particles, one of positively, 
and one of negatively charged particles. For spheres the inertia tensor J can be replaced by 
a single moment of inertia J. The total kinetic energy T is the sum of the kinetic energy, 
Ti, of the positive particles, and of the kinetic energy, T2, of the negative charges. 

We first perform the transformation above to the co-moving systems separately for the 
positive and the negative particles. This gives, 

T = ^M^Vl + \M2Vl + \hn\ + i J2O2 + T[ + T!,, (4) 

for the total kinetic energy, see Fig. [H In a second step we then introduce the co-moving 
system for the total system. We thus introduce, 

M = Ml + M2, J=Ji + J2, (5) 

/x = M1M2/M, / = Ji J2/J, (6) 

total mass and reduced mass as well as total moment of inertia and reduced moment of 
inertia. In terms of these one finds, 

V = (MiVi + M2V2)/M Vi = V + M2V/M (7) 

v = Vi-V2 V2 = V - Miv/M, (8) 



for the total center of mass velocity V, and relative velocity v, of the two spheres. Finally, 

n = (JiQi + J2^2)/J rii = n + Jiu/j (9) 

a; = fli - O2 0.2 = n- J2UJ/J, (10) 

gives the total average angular velocity fi, and the relative angular velocity cj, of the two 
oppositely charged spheres. In terms of these quantities we get the expression, 

T = ^MV^ + Ifiv^ + ^ JO^ + i/cj^ + r, (11) 

^ At ^ ^ 

for the total kinetic energy of our two-sphere system. The degrees of freedom in T' are 
assumed to be random and not to produce any net charge or current density. They will be 
ignored henceforth. 

III. LAGRANGIAN INCLUDING MAGNETIC INTERACTIONS 

Maxwell's equations and the equations of motion for the charged particles with the 



Is 



Lorentz force, can all be derived from a single Lagrangian via the variational principle |6|. 
The Lagrangian has three parts, particle, interaction, and field contributions. If radiation 
is neglected the field does not have independent degrees of freedom, but is determined by 
particle positions and velocities. Using the non-relativistic form for the kinetic energy one 
then gets, 

L = E {\n,v\ + 1^. . A(r,) - |0(r,)) , (12) 



where. 



and 



^) = Eir^> (13) 



\r — r 



A(r,t) = E^^^^^4%^. (14) 

Here the po.t,o,> a,.d ve.octy vectors of the part.e.es are ., a,.d .. resBecfve.y. ,„. a,.d , 



;heir masses and charges, and Cj = (r — rj)/|r — r^l (Darwin [7|, Jackson [8], Schwinger et al 



9|, Essen 10|, lUj)- The vector potential here is in the Coulomb gauge, and this essentially 
means that all velocity dependence of the interaction appears in the magnetic part, leaving 
the Coulomb interaction energy in its static form. 



When the expressions (fT3l) and iHM are inserted into equation (fT2|) one finds infinite 
contributions from self-interactions. When these are discarded, so that each particle only 
interacts with the field from the others, one obtains 

L = l^ 7,m,v, + ^ -^ '- -^ i '— - IJ -^, (15) 

where now r^j is the distance between particles i and ?' and Cj- is the unit vector pointing 



a 



from i to j. This is the so called Darwin Lagrangian [7|] for the system. We can write it, 

L = T + L^ag - $, (16) 

where, T, is kinetic energy and, $, the Coulomb electric interaction energy. The magnetic 
part can also be written 

^mag = E 1^^ ■ Mr^) = ^ / jXO " A{r) dV. (17) 

i 

Here it is important the the vector potential is divergence free (V ■ A = 0, Coulomb gauge). 
The Darwin Lagrangian thus includes both electric and magnetic interactions between the 
particles and is valid in when radiation can be neglected. 

IV. RELATIVE TRANSLATIONAL MOTION 

The Coulomb interaction, $, between two overlapping charged spheres is calculated in 
Appendix [XI The magnetic interaction between two charged spheres in relative translational 
motion is calculated in Appendix |B]for the case of small displacement of centers of the spheres 
{r <t^ R). Keeping only the quadratic term, one finds, 

L^^^ = 1^^2 i^^2 _ IQ! 2 (^g) 

Here r is the vector to center of the positive sphere from the center of the negative, so that 
r = V. The center of mass motion decouples, and we assume that the random motions 
decouple. This is then the relevant Lagrangian for the relative collective translation. It can 
be written. 



where, 



Lrel = \mv^ - \lCr\ (19) 



4Q2 / 4 A^re , 



and /C = Q^ / B?. If we assume a proton-electron plasma we get Mi = Nrrip, M2 = Nm, 
/i = Nmpin/{inp + m) ^ Nm, and Q^ = N'^e^ . On the right hand side we have introduced, 
Te = ;^, the classical electron radius. Note that when Nr^/R ^ 1 the effective mass Ai is 
entirely due to inductive inertia. 

Clearly the Lagrangian ( !T9l) corresponds to an oscillating system with angular frequency 



ojq = JK/M.. For this frequency we get explicitly, 

^0 = TfrwK- (21) 

If we introduce the dimensionless number, 

p ^ NrjR, (22) 

we see that for z/ ^ 1, one obtains, essentially, the Langmuir plasma frequency, u^ oc ne'^/m, 
see Eq. ([1]). If we reexpress the plasma frequency in terms of the classical electron radius 
[12], Eq. ( l2Ti) can be written in the form. 

Thus, when the number of particles is large enough so that z/ ^ 1, this gives. 

For this case the frequency turns out to depend on the size (radius) of the sphere, but not 
on the density. 



V. RELATIVE ROTATIONAL MOTION 

We now study pure rotational motion of the two charged spheres, about their coinciding 
centers of mass, but we include interaction, 

L, = -Jj{r)-A,{r)dV, (25) 

with a constant external magnetic field S = V x Ag. Here Ag = ^B x r is the vector 
potential of the external field. Starting from (ITTl) and ( IT6l) we find that, 

Lrot = 7:!'^ + -^mag + -^e, (26) 



is the relevant Lagrangian for collective relative rotation. 

The explicit calculations are sketched in Appendix O One finds that Lmag ~ ^^ and 
that this term therefore contributes to the effective moment of inertia, just as it contributed 
to the effective mass, Ai, in the translational case. The result can be written. 



where, 



L,ot = -Xa;2 + ^u: ■ B, (27) 

2 10c ^ ' 



2 Q'^R 2 . / 1 A^rc N 



We see that when u = Nr^/ R ^ 1 we can neglect the contribution from mass to the effective 
moment of inertia X. In this limit therefore, 

^ 2 Q^R 

^~— -^, 29) 

35 c2 ' ^ ' 

and there is essentially only inductive moment of inertia. We assume this below. 



VI. PLASMA ENERGY AND DIAMAGNETISM 

In order to investigate the equation of motion we put uj = ipe^ and B = i?(sin^e^ + 
cos^e^). The Lagrangian then becomes 

Lrot = ^XvJ^ + ^0 B cos e. (30) 

In general when, dL/dip = 0, one finds that, p^^ = {d/dt){dL/dip) = 0. In our case this 

gives, 

p^ = „ '^° = Tip H B cos 9 = const. (31) 

dip 10c 

If we assume that ip{t = 0) = when B{t = 0) = we find that the constant is zero: p<^ = 0. 
At all times we then find the relation, 

m = ^-^B{t)cose, (32) 

between angular velocity and the magnetic field. 

To get the energy from L{ip, ip) one calculates the Hamiltonian H = p^ip — L. For the 
Lj-ot of Eq. ( J30l) one finds. 



Let us consider two special cases of this phase space energy of the plasma. 

We first assume that the external field is zero {B = 0). The energy is then given by 
E = H = p'^/21. Here p^ is the angular momentum of relative rotation. For a given value 



of this angular momentum the energy is thus much smal 



er when z/ ^ 1 than otherwise. This 



reflects the fact, repeatedly stressed by the author [lO|, lul, ll3|, [lj|, that for given momenta 
the phase space energy of a plasma is lower when there is net current than in the absence 
of net current. 

Now consider instead the case p^ = 0. For simplicity we also assume ^ = 0. One then 
finds that 

or, equivalently, using f l29|) . that the energy as function of B is given by, 

^ ^ 80 40 V 3 ; Stt ^ ^ 

Note that here, B'^/Stt, is the energy density of the field B in our (gaussian) units. The energy 
is thus seen to grow quadratically with the applied magnetic field and our plasma spheres 
are strongly diamagnetic. Based on more detailed studies Cole [l5| has also concluded that 
plasmas are diamagnetic. In our model plasma diamagnetism is seen to be closely related to 
the diamagnetism of superconductors, as discussed by Essen [161]: the external field induces 
a current that screens the external field and reduces it inside. In the absence of resistance 
this screening current persists. 



VII. PLASMA RESISTIVITY 

Resistivity is completely neglected in the present model. It has been pointed out by 
Kulsrud [iTJ], in his book on astrophysical plasmas, that the negligible resistivity of such 
plasmas is in fact closely connected with magnetic induction. In the present treatment 
magnetic induction appears in the form of an inductive inertia that appears naturally as the 



main physical parameter in the present model. As early as 1933 Frenkel [18| suggested that 
superconductivity is due to inductive inertia. Frenkel also conjectured that inductive inertia 
can lower the energy and cause a phase transition. He did not, however, make his ideas 
quantitative. The present model system gives Frenkel's ideas some quantitative backing. 



We note that for a collective momentum p involving A^ particles the kinetic energy will be 
T = p'^ / {2mN) , when magnetic interaction is neglected, see Eqs. flT9l - l20|) . When the effect 
of magnetic interaction is included this becomes T + -Emag = p'^/{2Nm[l + 4i//5]) and we 
find that T + -Emag ^ T when u ^ 1, assuming that p remains constant. Collective modes 
are thus much more favorable thermodynamically when there is net current. 

Plasma resistivity is normally treated by studying the scattering of individual charged 
particles. Even with this type of treatment fast electrons become, so called, runaway elec- 



trons and experience no resistance [l9j. The present model indicates that resistivity can 
not be treated as resulting from the scattering of individual particles, since the collective 
motion of many charges leads to a large inductive non-local effect. All this points in the 
same direction, namely that plasmas need not be resistive, in agreement with our model 
treatment. 

If resistivity had to be included the translational oscillation would become a damped 
oscillation and any circulating current would eventually cease, thereby making the diamag- 
netic response temporary. A more immediate limitation of our model is probably the fact 
that a current would cause pinching and this would lead to instabilities that deform of the 
spherical shape. Lynden-Bell [22i] has studied the relativistically spinning charged sphere 
and finds that charge concentrates near the equator (as a result of pinching). 

VIII. CONCLUSIONS 

The model treated in this article is not particularly realistic. Instead it can be motivated 
as the simplest possible model within which one can study plasma phenomena associated 
with current, induction, and magnetic interaction energy, in a meaningful way. Hopefully it 



also has some novelty. In the literature one can find a 
oscillating plasma sphere (see e.g. Barnes and Nebel 



air amount of work on the radially 



3, 



Park et al. 



2l|), but not the 



modes treated here. A numerical study of a rotating convective plasma sphere, modelling a 



star, by Dobler et al. [23(1 shows how complicated more realistic models necessarily become. 
The two-sphere model studied here is therefore valuable as a device for gaining insight 
into some very basic plasma phenomena. As we have seen the most basic of these is the 
dominance of inductive inertia in the effective mass Ai of Eq. ( l20i) . and the effective moment 
of inertia X of Eq. (l28ll , when the number N of participating charged particles is large enough. 



One notes that v = Nvc/R ~ 118 for a typical laboratory plasma of density n = 10^° m~^ 
and of radius R = 1 cm, assuming that all particles contribute to the collective mode. Finally 
the model also indicates how this large inductive inertia influences the energy of the plasma 
and how a plasma responds diamagnetically to an external magnetic field. 



10 



Appendices 

Appendix A: Electrostatic interaction of t^vo overlapping charged spheres 

The electrostatic potential from a spherically symmetric charge density, g{r), is given by, 

0(r) = 4n (- r Q{ry ^ dr' + H ^(r>' dA . (Al) 

For a homogeneous charged sphere of total charge Q and radius R with g = Q/^^^, for 
r < R and zero for r > R, this gives, 

'i-S iorr<R 

(A2) 

— for r > R , 

and there is a pure Coulomb potential outside the sphere. 

Now consider a second, homogeneously charged, sphere of the same radius and of total 
charge q. If the distance r between the centers of the two spheres is greater than 2R the 
interaction energy is clearly given given by, 

aO 

$(r) = -i-^, for r > 2R. (A3) 

r 

We will now investigate what this interaction energy is when the distance is smaller so that 

the two spheres intersect. 

The interaction energy can be written 

^{r)= j Qq(t)dV (A4) 

where the volume integration is over the space occupied by the second sphere, where its 
charge density, Qq = g/^^y^, is different from zero. 

We do the integration using the slicing by concentric spheres depicted in Fig. [2l The 
volume element is then given by 

dV = t:-\R^ - ir - pf]dp. (A5) 

r 

When this is integrated between the limits p = r — R and p = r + R one should get the 

volume of the sphere and indeed and one easily checks that 

p=r+R 47]-_R3 

dV{p) = -^. (A6) 

11 



To do the integral (lA4p we have to spht the integration range at p = R since the function 
of flA2|) since it changes character at that radius. This gives 

rp=R f3Q Qp'^\ ^,, fp=^+RQ 



<^{r) = Qq 



dV + / — dV^ 

p=r-R\2R 2R^I Jp=R p 



(AT) 



The integrals are elementary and the final result, for r < 2R, is, 



$(r) = -^ (l92R^ - 80i?V2 + 30R^r^ - r^] . (A8) 

If we assume small displacements, r <t^ R, and put q = —Q, we find that, 

^^ 5R ^ 2i?3 ' ^ ^ 

so the potential is quadratic and there is there is a linear restoring radial force, F = —kr, 
with force constant, k = Q^/R^. 



Appendix B: Magnetic interaction of two charged spheres with relative translation 

Here we calculate Lmag, the middle term of Eq. flTSl) . for the two charged spheres with 
velocities Vi and V2 of Eqs. ([7]-[H])- This Darwin term in the Lagrangian gives explicitly, 

qi(l2 Y- Y- T^l ■ ^2 + (Vi ■ eik){V2 ■ Bik) , . 

^?? Fr^^^^ ^ '"" 

?2gi V- V- ^2 ■ Vi + (V2 ■ eki){yi ■ Bki) , . 

for the magnetic interaction. Here the indices i, j refer to the positive sphere and fc, / to the 
negative. The two terms f lB2|) and (]B3P representing interactions between the two oppositely 
charged spheres are clearly equal. If we denote the angle in fIBip between Vi and r^ — Vj 
by %, and similarly for ( JB4I) . we thus get 

, gi92 >^ Y- Vi • V2 + ( Vi ■ eik){V2 ■ Bik) 

12 



According to our assumptions the two double sums in (IBSP both represent the internal 
interaction between particles uniformly distributed within a sphere of radius R. They must 
thus be equal, and if we spht the Darwin Lagrangian into 

where Lgeif stands for f IBSP and Ljnt for ( ]B6|) . we find 

Without the squared cosines the double sum would simply be the Coulomb self interaction 
of a charged sphere, which can be taken from flA8l) with r = 0. Since the directions of 
Vi — Vj vary over the sphere it seems reasonable to estimate the effect of the squared cosine 
by replacing it with its spherical average 

1 /• o„.„ 1 



cos2 e = — cos' edn = -, (B9) 

An Jn 3 



where dQ = sin6'd6'd0. Also using. 



TV JV . r.nT2 

yV , , =^^, (BIO) 



we get 

/^2t/2 ^2t/2\ a a AT2 0/02 



.„.,4f,4f^i5^.|gl,vn^^, (BU) 



^ 4c2 4c2 / 35 i? 5i?c2 
where, Nqi = Q, Nq2 = —Q. We now introduce the transformation ©-([H]). If we put 

M2/M = a and Mi/M = /3 we can write it 

Vi = V + av, and V2 = V - f3v, (B12) 

where a + /3 = 1. This finally gives 



^self 



2^ l2V^ + 2{a-/3)V-v + {a^ + /3^y], (B13) 



5i?c2 
for the magnetic self interactions of the spheres. 

For the interaction part we put (]B12p into (]B6P and get 



r - Mi 



^, ^ 1 + cos^ _ ^^^2 ^ 1W^+ (B14) 






+ («-/3)E 



V- v + (V- eifc)(t;-eifc) 



i,fc 1^* ^*:l 
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(B15) 



We now split the vector v into a part parallel to V and a part perpendicular to V, to get 
V = vy + v±. The double sum in flB15|) is then 

.^V -v + iV ■ eik){v ■ eik) ^-^ 1 + cos"^ 9 ik , ., ,^ cos 6*^^ cos 6* -fc 

Z^ \ — :: \ = ^^v2^—, — r- + Vv±2_^ — ■ — ■ — . (B16) 

In the above expressions, 9ik, 0[^, and 9'/i^, are the angles between rj — Vk and V, v, and v±, 
respectively. The double sums are now over points distributed in displaced spheres. If we 
assume that the displacement r is small compared to the radius R we can again approximate 
the effect of the cosines by spherical averaging flB9p . We then find that 



l + cos2^= -, and, cos^cos^' = 0, (B17) 



where 6 and 6' represent angles with perpendicular directions. 
According to our assumptions we can now use Eq. ( 1A9I) and put 



i-,w^r^-w^- (BIS) 

Finally then we get (note that Vvy = V ■ v) 

L-., = -^ (^ - ^^) W + 2(« -m-v- 2aPv\ (B19) 

for the magnetic interaction of the two spheres. 

If we now add the self interaction, Eq. ( JB13I) to the result ( 1B19I) we end up with. 



'""^ ~ 10 i? C2 + 3 i? i?2 c2 ' ^ ^> 

for the total magnetic interaction part of the Lagrangian. Here we have used, a + /3 = 1. 
Note also that we have assumed that r <^ R and that certain angular dependencies have been 
treated approximately. If we also assume that v"^ <^ c^ the second term here is negligible 
compared to the first. 

Appendix C: Magnetic interaction of t^vo charged spheres ^vith relative rotation 

If we put Qi = Q/^^Y^ and ^2 = —^1, for r < R, and zero outside, the two spheres have 
current densities, 

j,(r) = ^,axr, z = l,2. (CI) 

14 



Divergence free vector potentials, that match smoothly with dipole field vector potentials 
outside the spheres, are then given by [24|, 

A.(r) = |(^i?^-r^)j.(r). (C2) 

Using this, it is elementary to show that, 

i^mag = ^ / J(^) ■ Mr) ^y = Yc /(-^'i ■ ^1 + -^'2 • M + 2Ji ■ A2) dV^, (C3) 

is given by, 

L^., = |^-^ (C4) 

where, uj = Q,i — 0,2-, is the relative angular velocity. 

The response to a constant external field B with vector potential A^ = \B x r, is given 
by, 

Le = -j3{r) ■ A,{r)dV = -J{3i ■ A^ + j^ ■ A,) dV. (C5) 

One finds, 

L, = 9^B ■ uj, (C6) 

10c ' ^ ' 

from straightforward calculations. 
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Figure captions 

FIG. 1: Some notation for the two-sphere plasma model treated in this article. Note the separation 
of the spheres is exaggerated for clarity. It is assumed small in the relative translation case and 
zero in the relative rotation case. 

FIG. 2: Integration over the upper sphere, centered on q, is performed by means of volume elements 
dV = Sdp that consist of the space between two concentric sphere segments centered on Q with 
radii p and p + dp. The area S of such a segment is given by S" = 2Tr ph where h = p{l — cos a). 
The cosine theorem applied to the triangle QqP gives B? = r''^ + p^ — 2rpcosa. This gives cos a 
and insertion into d^ gives dV = {iTp/r)[E? — (r — p)'^]dp. 
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